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The dynamical behavior of propagating structures, determined from a Karhunen-Loeve decomposi- 
tion, in turbulent pipe flow undergoing reverse transition to laminar flow is investigated. The tur- 
bulent flow data is generated by a direct numerical simulation started at a fully turbulent Reynolds 
number of Re r = 150, which is slowly decreased until Re T = 95. At this low Reynolds number the 
high frequency modes decay first, leaving only the decaying streamwise vortices. The flow under- 
goes a chugging phenomena, where it begins to relaminarize and the mean velocity increases. The 
remaining propagating modes then destabilize the streamwise vortices, rebuild the energy spectra, 
and eventually the flow regains its turbulent state. Our results capture three chugging cycles be- 
fore the flow completely relaminarizes. The high frequency modes present in the outer layer decay 
first, establishing the importance of the outer region in the self-sustaining mechanism of wall bound 
turbulence. 



I. INTRODUCTION 

In part I the effect of drag reduction on propagat- 
ing turbulent flow structures by spanwise wall oscillation 
was studied. 1 A second instance where drag reduction is 
seen is in a relaminarizing flow. As the turbulence dies, 
so does the Reynolds stress generation, and thus, for a 
constant pressure gradient driven flow, the flow rate in- 
creases. As we will show, the flow does not immedi- 
ately relaminarize, but instead goes through a series of 
chugging motions. In these chugging motions, the flow 
loses its turbulent inert ial range, losing the high frequen- 
cies first. Before the flow has completely relaminarized, 
certain key propagating waves interact with the decay- 
ing streamwise vortices, recreating the cascading energy 
scales that populate the inertial subrange. In this part, 
we examine the dynamics found in relaminarization to 
understand how a flow remains turbulent to better eluci- 
date the mechanism behind the self-sustaining nature of 
turbulence. 

Previous work has focused either on relaminar- 
ization from favorable pressure gradients j2*£d strong 
accelerations 5 , or examining relaminarization as a 
testbed to understand the decay rate of structures^ Out- 
side the field of wall turbulence, there has been related 
work in studying relaminarization and bifurcations in 
spherical Couette flow. 7 The only work found that relates 
to drag reduction is the examination of linear feedback 
control in a turbulent channel flow that achieved total 
relaminarization. 8 Nevertheless, it is of interest to exam- 
ine the field of transition where recent work has added 
to our knowledge of pipe flow structures and their inter- 
actions. 

Both Kerswell 9 and Faisst and Eckhardt 10 have found 
traveling wave solutions to the Navier- Stokes equations 
through continuation methods. They identified struc- 
tures for rotationally symmetric solutions, which is con- 
firmed here for Re r = 95 and in previous workii for 
Re r = 150 through a Karhunen-Loeve (KL) decomposi- 
tion of a direct numerical simulation of turbulent pipe 
flow. Moreover, Kerswell found that the threefold rota- 



tion, or azimuthal wavenumber n = 3 in our notation, is 
the largest contributor near the critical bifurcation point 
associated with the laminar to turbulent transition. This 
also corresponds to our findings, as we will show that 
when observing the chugging phenomena, the n = 3 trav- 
eling wave is often the most energetic at the point which 
the flow reasserts itself as turbulent. 

Also of note is the minimal channel work by Webber, 
Handler, and Sirovich; we apply their results to further 
the understanding of this chugging phenomenal They 
indicate that the nonlinear terms in the Navier-Stokes 
equations lead to triad interactions of the KL modes 
which is responsible for the transfer of energy between 
modes. This occurs whenever wavenumbers of three 
modes (m, n, g), (ra', n', q') : and {m" n" , q") sum to zero, 
shown in equations [1] and [2] below, 



n + ri + n" = 
m + m! + m" = 



(i) 

(2) 



where n is the azimuthal wavenumber and m is the 
streamwise wavenumber obtained from the Fourier rep- 
resentation of the flow. 



II. NUMERICAL METHOD 

Details on the numerical method for generating the 
direct numerical simulation (DNS) flow fields and on the 
Karhunen-Loeve method can be found in part I of this 
paper. 1 

Over a time of 12,000 t + the Reynolds number was 
slowly reduced from Re r = 150 to a value of Re r = 95, 
and the DNS continued for another 5000t + to eliminate 
any transient effects, as seen in Figure [TJ Data was 
then collected for 10,000t + , which included three dis- 
tinct chugs before the flow completely relaminarized, as 
seen in Figure [2j 

The grid resolution was kept the same as for the orig- 
inal Re r = 150 case. Thus the grid is effectively further 
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FIG. 1: Time history of the mean flow rate relaminarization 
showing the initialization from Re T = 150 to Re r = 95. The 
window shows where data was collected. 
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FIG. 2: Time history of the mean flow rate relaminarization 
for Re T = 95. Three chugging cycles are seen near t + = 
2000, 4000 and 6000 with the final relaminarization starting 
near t + = 8000. 

refined to Ar+ « 0.49 and (RAQ)+ « 3.1 near the wall 
and A + ^2.0 near the centerline with a constant stream- 
wise resolution of Az + = 4.0, where r, 0, and z is the 
radial, azimuthal, and streamwise direction, respectively, 
and R is the radius of the pipe. 



III. RESULTS 

The profile of the mean flow with respect to wall units 
is seen in Figure [3j At this low Reynolds number, the 
profile does not conform to the log layer, yet near the 
wall it still exhibits a linear sublayer. Also, the flow does 
not conform to the laminar parabolic profile, indicating 
that it is indeed turbulent. 

The root-mean-square (rms) velocity profiles also show 
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FIG. 3: Mean flow profile versus y + — (1 — r)Re r including 
the theoretical turbulent profile (dashed) with the sublayer 
(u + — the log layer (u + — log(^/ + )/0.41 + 5.5), and the 
analytical parabolic laminar solution (dash-dot). The mean 
flow profile follows the law of the wall, yet deviates drastically 
from the log layer as expected. 
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FIG. 4: Root-mean-squared velocity profiles for Re T = 95 
(solid) and Re T = 150 (dashed) versus wall units y + . The 
inflections near y + — 90 are effects of the laminar chugging, 
because the rms velocities are averaged over all time steps. 

a turbulent trend as seen in Figured! although in compar- 
ison to the Re r = 150 flow, the radial and azimuthal fluc- 
tuations are about half as strong, and the peak stream- 
wise rms velocity is shifted further away from the wall. 
Also noteworthy is the strong fluctuations near the cen- 
ter of the pipe at y + ~ 90. This is the first indication 
that the dynamics near the center of the pipe differ from 
what is expected for the fully turbulent case. 

The Reynolds stress profile also differs from the Re r = 
150 case and is shown in Figure The Reynolds stress 
for the Re r = 95 case has roughly half the magni- 
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FIG. 5: Reynolds stress profiles for Re T = 95 (solid) and 
Re T = 150 (dashed) versus wall units y + . As in the rms 
velocity profiles, the inflections near y + = 90 are effects of 
the laminar chugging, showing that relaminarization begins 
at the center of the pipe. 



FIG. 6: Comparison of the running total energy retained 
in the KL expansion for Re T = 95 (solid) and Re T = 150 
(dashed). The 90% crossover point contains 2453 and 66 
modes respectively, showing a drastic reduction in the tur- 
bulent attractor. 



tude throughout the radial profile, although the peak 
Reynolds stress is also found at the same location of 
y + = 31. The largest deviation from the expected profile 
occurs near the centerline beyond y + ~ 60, where the 
Reynolds stress begins to fluctuate. In particular, be- 
tween y + = 88 and 92, the Reynolds stress is negative, 
which is physically interpreted as turbulence damping. 
Now, in addition to the rms deviation near the center- 
line, this Reynolds stress fluctuation indicates that the 
relaminarization process begins at the center of the pipe 
and goes towards the wall. 

Turning from statistics to the KL decomposition, we 
find that the chaotic attractor is reduced in size, as ex- 
pected with a reduction in Re r , from Dkl — 2453 to 
Dkl = 66, shown in Figure [6j This dimension is simi- 
lar to that found in Part I where the oscillated pipe was 
barely turbulent with a dimension of Dkl = 102, and 
any stronger oscillation would have resulted in relami- 
narization. 

In observing the energy content of the modes in Table 
[T| and the most increased and decreased in Table [Til we 
find that the shear modes increase drastically, which is a 
result of the chugging motion and large mean flow rate 
fluctuations. Also of note is the increase in strength of 
the n = 3, 4 and 5 streamwise vortices (m = 0) and their 
associated wall traveling waves (m = 1). The increase in 
the (0, 1, 1) and (1, 0, 1) modes, since they are not found 
in the work by Kerswell 9 and Faisst and Eckhardt, 10 
could be the catalysts of energy triads in equations [1] and 
[2] between the n = 3,4 and 5 rolls and the wall traveling 
waves. For example, the (1, 3, 1) and (1, 4, 1) waves inter- 
act through the (0,1,1) catalyst, and the (1,3,1) wave 
and the (0, 3, 1) roll interact through the (1, 0, 1) catalyst. 
The modes that decreased the most in energy are more of 
a result of the lower Reynolds number, as the high modes 
n = 8, 9 and 10 in the Re r = 150 case were an important 
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FIG. 7: Comparison of normal speed locus for the Re T = 95 
(•) and Re r = 150 (+) cases. The solid lines represent a circle 
of diameter 8.41 and 8.64 respectively that intersect at the 
origin. 



basis for streamwise vortices that were 94 — 117 wall units 
apart. However, at Re r = 95, these represent spacings 
of 60 — 75 wall units, which is too small, thus explaining 
their drop in energy. 

Looking at the normal speed locus in Figure El the 
modes for a wave packet similar to that of Re r = 150, 
but with a slightly faster advection speed of 8.64 versus 
8.41. This shows a Reynolds number dependence on the 
advection speed, as expected, since the advection speed 
was shown to scale with the mean flow rate when the flow 
rate was increased in part I with spanwise wall oscillation. 

In visualizing the most energetic modes, the same 
structures as those found in the Re r = 150 case are 
present. Figures [8] - [13] show the coherent vorticity for 
the most energetic modes for each of the subclasses dis- 
cussed in part I. Since the (0, 0, 1) mode has no coherent 
vorticity, the velocity is shown. Here we use "coherent 
vorticity" to refer to the imaginary eigenvalues of the ve- 
locity gradient tensor as defined in Chong et al. (also 
know as "degree of swirl"). 13 Although slight differences 
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TABLE I: Comparison of first 25 eigenvalues, m is the stream- 
wise wavenumber, n is the spanwise wavenumber, and q is the 
eigenvalue quantum number. 



Re T = 


150 


Re T = 95 


Index m n q Ener£ 


;y % Total 


m n q Energy % Total 



x 


o 


6 


1 


1.61 


2.42% 


o 


o 


x 


114 


55 76% 


2 


n 
u 


u 




1.48 


2.22% 


n 

w 


o 


x 


y.uo 


4 79% 


Q 
o 


n 
u 


Q 
o 




1 4^ 


9 1 7% 

£i.X. i /0 


n 


4 


x 


o.oy 


4 1 1 % 


4 


o 


4 




1.29 


1.93% 


o 




x 


7.79 


3.82% 


5 


o 


2 


I 


1.26 


1.88% 


o 


2 


x 


6.30 


3.09% 


6 


X 


5 


I 


0.936 


1 40% 

± .'-tKJ /{J 


o 


o 


2 


5.98 


2.93% 


7 


X 


6 


I 


0.917 


1 S7% 


o 


5 


x 


2.79 


1 S7% 


g 


X 


3 


I 


0.902 


1.35% 


I 


4 


x 


2.18 


1 07% 


q 


X 


4 


I 


0.822 


1 91% 


I 


o 


x 


1 

J- .ou 


QD% 

U.C7U /0 


i n 


n 
u 


I 


I 


u.ouu 


1 9f)% 

-L .ZiU /0 


n 


n 


o 


J- .oo 


QD% 

U.C7U /0 


1 1 
11 


1 
1 


( 


l 


U. / Do 


1 1/10/ 

1.14/0 


U 


O 


1 


1.00 


U. M /O 


12 


1 


2 


1 


0.683 


1.02% 


1 


5 


1 


1.45 


0.71% 


13 





7 


1 


0.646 


0.97% 


1 


2 


1 


1.42 


0.70% 


14 


2 


4 


1 


0.618 


0.92% 





1 


2 


1.34 


0.66% 


15 





8 


1 


0.601 


0.90% 


1 





1 


1.30 


0.64% 


16 


2 


5 


1 


0.580 


0.87% 


1 


6 


1 


0.884 


0.43% 


17 


1 


1 


1 


0.567 


0.85% 


1 


1 


1 


0.815 


0.40% 


18 


2 


7 


1 


0.524 


0.78% 





7 


1 


0.727 


0.36% 


19 


1 


8 


1 


0.483 


0.72% 


2 


4 


1 


0.716 


0.35% 


20 


2 


6 


1 


0.476 


0.71% 


2 


3 


1 


0.665 


0.33% 


21 


2 


3 


1 


0.454 


0.68% 


2 


5 


1 


0.589 


0.29% 


22 


2 


2 


1 


0.421 


0.63% 


2 


2 


1 


0.544 


0.27% 


23 


2 


8 


1 


0.375 


0.56% 





2 


2 


0.521 


0.26% 


24 


1 


9 


1 


0.358 


0.54% 


2 


6 


1 


0.487 


0.24% 


25 


3 


4 


1 


0.354 


0.53% 


1 


1 


2 


0.479 


0.23% 



TABLE II: Ranking of eigenfunctions by energy change be- 
tween the Re T = 150 and Re T = 95 cases, m is the stream- 
wise wavenumber, n is the spanwise wavenumber, and q is 
the eigenvalue quantum number. The shear modes (0, 0, q) 
are among the most increased, as are the n = 3, 4, and 5 rolls 
(m — 0) and wall modes (m = 1). 
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are visible due to the lower Reynolds number, the same 
trends and characteristics are found for the KL modes as 
those found in the Re r = 150 case. 
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FIG. 8: The (0,0,1) shear mode with contours of stream- 
wise velocity and vectors of cross-stream velocities. Left: (a) 
Re T = 95. Right: (b) Re T = 150. 
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FIG. 9: The (0,3,1) roll mode with contours of coherent vor- 
ticity. Left: (a) Re T = 95. Right: (b) Re T = 150. 
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FIG. 10: The (1,4,1) wall mode with contours of coherent 
vorticity. Left: (a) Re T = 95. Right: (b) Re T = 150. 
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FIG. 11: The (3,2,1) lift mode with contours of coherent vor- 
ticity. Left: (a) Re T = 95. Right: (b) Re T = 150. 
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FIG. 12: The (1,1,1) asymmetric mode with contours of co- 
herent vorticity. Left: (a) Re T = 95. Right: (b) Re T = 150. 



IV. DYNAMICS 

Recreating the time history of the KL modes reveals 
the interaction between the shear modes, roll modes, and 
propagating waves. As shown in Figure [HJ the chugging 
phenomena happens when the propagating modes drop in 
energy. This happens at times t + « 1600, t + « 3500, and 
t + « 5500, and starts the chugging phenomena. At t + « 
8000, the propagating modes drop too far in energy, and 
the flow relaminarizes. The chugging cycle ends again 
when the propagating waves spike at t + « 2200, t + « 
4100, and t+ « 6200. 

Examining the propagating waves based upon their 




FIG. 13: The (1,0,1) ring mode (r-z plane cross-section) with 
contours of coherent vorticity. top: (a) Re T = 95. Bottom: 
(b) Re r 150. 




FIG. 14: Time history of energy of shear modes (blue), roll 
modes (black), and propagating waves (red). 
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FIG. 15: Time history of energy of propagating subclasses. 
Wall (solid black), lift (red), asymmetric (dashed), and 
rings (blue). 




FIG. 16: Ratio of lift mode energy to wall mode energy as 
a function of time. As the lift mode energy decreases with 
respect to the wall mode energy, relaminarization begins. 



subclass, shown in Figure[15l we find that the asymmetric 
and ring modes stay about the same energy as each other 
throughout the chugging cycles, while the wall modes are 
about a factor of 4 times more energetic. The lift modes, 
on the other hand, vary greatly in energy. When the 
lift modes decay from high to low energy, this coincides 
with the start of a chugging cycle, and ends when the lift 
modes regain the high energy state. This is emphasized 
in Figure [16] where the ratio of the lift mode total energy 
to the wall mode total energy is shown. In Figure \T7\ 
only one chug is shown, emphasizing the importance of 
the lift to wall ratio, and also the energy transfer between 
classes as energy flows from the rolls to the wall modes 
(through the ring modes), and then from the wall modes 
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FIG. 17: Time history of energy of subclasses for a single chug. 
From the top: the mean velocity, roll mode energy, ring mode 
energy, wall mode energy, asymmetric mode energy, lift mode 
energy, and the lift to wall mode energy ratio. The start of the 
chug (blue dashed line) is when the lift to wall mode energy 
ratio drops too low. The end of the chug (red dashed line) is 
when the ratio spikes and recovers. The phase lag in energy 
from the roll to the wall modes (through the ring modes) and 
similarly from the wall to the lift modes (primarily through the 
asymmetric modes) can be seen. 
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FIG. 18: Average energy spectra of the propagating modes. 
Wall (solid), lift (dashed), asymmetric (dots), and rings (dash- 
dot). 
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FIG. 19: Time history of the energy spectra of the propagating 
modes. t + = and t + = 2392 show the proper established 
energy spectra at the beginning of the simulation, and well 
after the chug. t + = 1592 is the spectra when the chug begins. 
t + = 2112, 2552 and 2232 show the spectra regaining strength 
with time, finishing the chug cycle. Like the dynamics of the 
lift to wall energy ratio, the chugging spectra reinforces the 
dependence of the self-sustaining mechanism of wall turbulence 
on the high frequencies found in the inert ial range, represented 
by the lift modes. 



to the lift modes (primarily through the asymmetric 
modes). 

Observing the energy spectra throughout the chug, we 
first revisit the total energy distribution of the propa- 
gating waves as averaged over the entire flow, seen in 



Figure [TBI Again, like the Re r = 150 spectrum, the lift 
modes are more energetic than the wall modes for high 
wavenumber s. Thus, because the high wavenumber s of 
the spectra dies off by two orders of magnitude at the 
start of the chug cycle, seen in Figure [HE it reinforces 
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the importance of the lift modes in maintaining the tur- 
bulent flow. 

As noted in the Re r = 150 case, the lift modes are re- 
sponsible for the majority of turbulence near the center 
of the pipe, as the wall modes stay near the wall, even 
for high quantum number. Thus, the importance of the 
outer region in the self sustaining mechanism of turbu- 
lence is reinforced. If the lift modes do not receive enough 
energy, cascaded through the wall modes, the relaminar- 
ization process begins near the center of the pipe. This 
is confirmed by the presence of the large drop in the rms 
velocities and Reynolds stress profiles near the center of 
the pipe. If the process cannot be halted in time by the 
transfer of energy from the wall modes to the lift modes, 
typically through the (1,3,1) or (1,4,1) modes, the flow 
will completely relaminarize. Adding in the findings of 
triad interactions by Webber et al.^ the flow of energy 
is shown in Figure [20j The shear to roll interactions are 
catalyzed by the rolls themselves, the roll to wall inter- 
actions are catalyzed mostly by the ring modes, and the 
wall to lift modes are catalyzed by both the ring and the 
asymmetric modes, and other lift modes. In the relami- 
narization process, it is this final leg that fails, breaking 
the mechanism, and starting relaminarization from the 
center of the pipe. 



CONCLUSIONS 



It is apparent through the examination of the spanwise 
wall oscillated case in Part I and the relaminarization 
case in Part II that if any leg of the energy cycle in a 
turbulent flow is disrupted, the resulting imbalance can 
lead to the start of a relaminarization process, and even 
complete relaminarization. In Part I, we describe a model 
where the propagating wall modes were pushed to higher 
advection speeds, reducing their effective lifespan. Thus 
they do not have enough time to take energy from the roll 
modes, breaking the third leg of the mechanism. Here in 
Part II, with lower pressure gradient, there is not enough 
energy to properly maintain the lift modes, and the last 
leg of the process is broken, starting the relaminarization 
process. 

Thus, in conclusion, we find that while the wall modes 
and near wall interactions are responsible for the genera- 
tion of turbulence from the pressure gradient, the turbu- 
lence in the outer region is necessary to maintain the 
proper inertial range in the energy spectra, and that 
without it, the relaminarization process begins. 
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